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Abstract  -  In  this  paper  we  investigate  the  effect  of  an  inelastic  interface  layer  on  the 
mechanical  behavior  of  a  transversely  loaded  fiber-reinforced  composite.  A  simple  linearly 
viscoelastic  model  is  used  to  characterize  the  stiffness  and  viscosity  of  the  interface  separating 
the  fiber  and  matrix  phases.  The  mechanical  response  is  obtained  using  the  finite  element 
method  and  calculations  are  carried  out  for  a  unit  cell  in  a  periodic  array  of  hexagonally  packed 
fibers.  An  approximate  representation  of  the  time  dependent  macroscopic  behavior  of  the 
composite  is  derived  analytically  and  compared  with  the  numerical  results.  From  a 
micromechanical  perspective,  tha  mf'.c'nce  of  interfacial  stress  relaxation  on  the  stress  fields 
in  the  matrix  material  contiguous  to  the  interface  is  examined. 


1.  INTRODUCTION 

Since  the  advent  of  the  modem  fiber-reinforced  composite,  the  determination  of  the 
mechanical  properties  of  these  materials  has  become  of  significant  practical  importance. 
Unlike  the  axial  strength  and  stiffness  properties  which  are  primarily  governed  by  the  axial 
properties  of  the  fiber,  the  behavior  of  the  fiber-reinforced  composite  in  the  transverse 
direction  is  dominated  by  a  relatively  low  stiffness  matrix  material  and  the  fiber/matrix 
interface.  This  may  place  severe  limitations  on  the  overall  performance  of  the  composite 
and  thus  it  is  desirable  to  accurately  characterize  the  transvers  operties. 

In  most  analytical  and  numerical  work,  investigators  have  assumed  a  perfect  bond 
between  the  fibers  and  the  matrix  material  which  is  modeled  by  continuity  of  interfacial 
tractions  and  displacements.  In  reality,  however,  a  more  complex  state  exists  between  thd' 
fiber  and  matrix  constituents,  and  the  assumption  of  perfect  bonding  may  not  be  suitable  in 
the  presence  of  a  thin  interfacial  zone  which  connects  the  two  phases  (e.g.  fiber  coating  or 
intermolecular  bonding).  In  this  analysis  it  is  assumed  that  the  interface  layer  is 
infinitesimally  thin  and  supports  a  traction  field  having  both  normal  and  tangential 
components.  Continuity  of  tractions  is  assumed  across  the  interface,  however,  fiber  and 
matrix  phase  separation  is  simulated  by  allowing  for  displacement  discontinuities  across  the 
interface.  Such  a  model,  assuming  a  linear  relationship  between  the  displacement 
difference  across  the  interface  and  the  conjugate  tractions,  is  employed  by  Aboudi  (1987), 
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Steif  and  Hoysan  (1987),  Achenbach  and  Zhu  (1989a, 89b)  and  Hashin  (1989). 
Needleman  (1987)  utilizes  this  model  to  simulate  the  bond  between  rigid  spherical 
inclusions  embedded  in  an  isotropically  hardening  elastic -viscoplastic  matrix  and  considers 
a  more  general  interfacial  constitutive  relation. 

In  implementing  the  above  mentioned  interface  model,  care  must  be  taken  in  order  to 
avoid  an  unrealistic  interpenetration  of  the  matrix  and  fiber  phases  which  can  occur  in  local 
regions  of  compression.  As  will  be  seen  later,  the  imposition  of  such  a  constraint 
significantly  influences  the  transverse  mechanical  behavior  of  the  composite  with  a 
relatively  low  stiffness  interfacial  zone. 

In  the  present  paper,  the  role  of  an  inelastic  interface  on  the  transverse  properties  of 
the  fiber-reinforced  composite  is  considered  from  both  a  macroscopic  and  a  microscopic 
perspective.  The  choice  of  a  relatively  simple  linearly  viscoelastic  interface  model  allows 
us  to  derive  an  approximate  representation  for  the  transverse  relaxation  moduli  of  the 
composite  with  which  to  compare  our  numerical  results.  We  note  that  this  analytical  model 
does  not  incorporate  the  aforementioned  impenetrability  constraint,  and  of  primary  interest 
is  how  well  it  compares  with  the  numerical  results  over  a  lange  of  interfacial  stiffness 
parameters.  In  addition,  the  influence  of  interfacial  stress  relaxation  on  the  stress  fields  in 
the  matrix  material  contiguous,  to  the  interface  is  examined  from  a  micromechanical 
viewpoint. 

In  the  following  section,  the  model  which  is  chosen  to  represent  the  behavior  of  a 
unidirectional,  fiber-reinforced  composite  is  described  and  the  boundary  value  problem  is 
formulated.  Discussed  in  section  3  are  implementational  details  of  the  finite  element 
method,  the  numerical  procedure  employed  in  this  analysis.  In  section  4,  the  analytic 
model  which  approximates  the  transverse  relaxation  moduli  of  the  fiber-reinforced 
composite  is  described  and  compared  with  the  numerical  results.  In  section  5,  a  qualitative 
discussion  is  presented  regarding  the  effect  of  interfacial  stress  relaxation  on  the  resulting 
stress  fields  in  the  matrix  material  near  the  fiber/matrix  interface.  Finally,  some  concluding 
remarks  are  stated  in  section  6. 


2.  FORMULATION  OF  THE  BOUNDARY  VALUE  PROBLEM 

A  cross-sectional  view  of  the  model  employed  in  this  analysis  is  illustrated  in  figure  1. 
It  is  assumed  that  the  fibers,  all  of  equal  radius,  a,  are  periodically  spaced  in  a  regular 
hexagonal  array  and  are  '“■"bedded  in  an  infinite  matrix.  Two  loading  directions  are 
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considered  in  this  analysis,  the  closest  packing  direction  (CPD)  and  the  mid-closest  packing 
direction  (Mid-CPD).  The  Mid-CPD  loading  direction  which  bisects  the  angle  formed  by 
two  closest  packing  directions  is  illustrated  in  figure  1. 

Through  arguments  of  symmetry,  it  is  only  necessary  to  analyze  the  rectangular  region 
outlined  in  figure  1  and  shown  in  detail  in  figure  2.  Neglecting  all  rigid  body  motion,  the 
point  O,  located  at  the  origin  of  the  Cartesian  coordinate  system  shown  in  the  figure,  is 
considered  fixed  throughout  the  analysis.  By  imposing  the  appropriate  boundary' 
conditions  on  one  half  of  rbis  rperan«nlar  region,  i.e.  the  trapezoid  AEEF  consisting  of  one 
quarter  of  a  regular  hexagon  with  sides  of  length  b,  the  state  of  stress  and  strain  for  the 
entire  model  may  be  completely  characterized. 

2.1.  Boundary  Conditions 

The  relevant  boundary  conditions  for  the  case  of  CPD  loading  are  given  below.  The 
corresponding  expressions  for  the  case  of  Mid-CPD  loading  are  obtained  similarly,  see  e.g. 
Achenbach  and  Zhu  (1989b).  Referring  again  to  figure  2,  the  loading  direction  is  parallel 
to  the  X2  axis  of  the  Cartesian  reference  frame  centered  at  O.  Relative  to  this  frame,  the 
boundary  conditions  along  the  external  boundaries  AB,  EF,  and  AF  respectively  are 
expressed  as  follows: 
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where  u/  and  u2  are  the  displacement  components  in  the  X/  and  j <2  directions  and  A2  is  the 
magnitude  of  the  prescribed  displacement  in  the  x2  direction  along  AB  and  EF.  The 
quantity  Aj  is  the  magnitude  of  the  unknown  displacement  in  the  X/  direction  along  AF  to 
be  determined  as  part  of  the  numerical  solution.  Along  BE,  the  following  displacement 
condition  must  hold: 


ut(-xi,-x2 )  =  u  [{x  1  yc2),  u2(-xi-x2)  =  u2(xhx2).  (4) 

One  additional  relevant  condition  for  this  case  of  loading  is  obtained  through  equilibrium 
considerations  in  the  x2  direction  and  is  stated  as 
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where  T?  is  the  traction  component  in  the  X2  direction  on  the  external  boundaries  AB  and 
EF.  This  condition  allows  for  the  numerical  determination  of  the  remote  applied  stress  <x». 

2.2.  Constitutive  Relations 

Now1  that  the  relevant  conditions  that  must  be  satisfied  on  the  external  boundaries  of 
the  trapezoidal  region  ABEF  have  been  given,  we  turn  our  Mention  to  the  individual 
material  phases  comprising  the  interior  of  the  trapezoidal  region.  In  the  present  analysis,  a 
graphite/epoxy  composite  material  system  is  considered,  and  a  description  of  the 
constitutive  law  which  governs  the  behavior  of  each  phase  follows. 

Epoxy  Matrix  and  Graphite  Fiber  Phases.  It  is  assumed  that  the  matrix  is  isetmpic 
and  linearly  elastic.  The  fiber  is  taken  to  be  linearly  elastic  and  transversely  isotropic.  The 
elastic  constants  employed  in  this  analysis  were  obtained  by  Kriz  and  Stinchcomb  ( 1979; 
and  are  given  in  table  1.  For  the  case  of  plane  strain,  the  sD'ess-strain  relations  for  the 
matrix  phase  can  be  written  as 

t '33  =  £JJ  -  £23  -  0  (6) 
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where  the  parameters  X  and  p  are  the  Lame'  constants  and  the  Greek  indices  a,  [3,  and  y 
range  over  1  and  2.  For  the  transversely  isotropic  fiber  phase,  the  elastic  stress-strain 
relations  for  the  case  of  plane  strain  are  given  by  Hashin  (1979).  The  non-zero 
components  are  written  as 


Ojj  =  l  (£,i  +  £22) 


(7) 
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where  the  constant  /  is  related  to  the  axial  Poisson's  ratio, v^,  and  the  transverse  bulk 
modulus,  Kj,  by 


/  =  2KTvA 


(8) 


In  addition,  the  transverse  hulk  modulus,  Kj ,  is  related  to  the  engineering  constants  by  the 
relation 
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where  Ej  and  Ea  are  the  transverse  and  axial  Young's  moduli  of  the  fiber  and  Gy  is  the 
transverse  shear  modulus. 


Interface  Model.  Both  a  linearly  elastic  and  a  linearly  viscoelastic  constitutive  relation 
are  considered  for  the  interfacial  zone.  For  the  linearly  elastic  interface,  it  is  assumed  that 
the  normal  traction  between  the  fiber  and  matrix  phases  is  proportional  to  the  corresponding 
normal  relative  dtspk  cement  across  the  interface.  Similarly,  the  tangential  traction  is  taken 
to  be  proportional  to  the  tangential  relative  displacement  across  the  interface.  Thus, 

Tn  =  kn[un],  (10) 


Tt  =  ktlu,], 


where 


un  =  Uintn  and  Tn  =  ( Cijmnfn  (11) 

ut  =  u  -  un  and  Tt-T-Tn 

are  the  normal  and  tangential  displacement  and  traction  vectors  respectively,  [  .  ]/  denotes 
the  jump  in  the  relative  quantity  across  the  interface,  kn  and  kt  are  normal  and  tangential 
stiffness  parameters  arbitrarily  taken  to  be  equal  in  this  analysis,  and  T  is  the  traction 
vector  (7j  =  Here  Ofy  =  oJt  is  the  Cauchy  stress  tensor.  Since  we  are  considering 

the  plane  strain  problem,  the  indices  i  and  j  in  (1 1)  take  on  the  values  1  and  2  referring  to 
the  Cartesian  coordinate  system  shown  in  Figure  2.  A  positive  normal  relative  displacement 
[«„]/  denotes  normal  separation  between  the  fiber  and  matrix  phases.  However,  we 
assume  that  a  negative  normal  relative  displacement  [w„]/  allows  for  a  physically  unrealistic 
interpenetration  of  the  matrix  phase  into  the  fiber  phase  and,  thus,  we  enforce  the 
impenetrability  constraint 
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(12) 


For  the  linearly  viscoelastic  interface,  the  time  dependent  response  of  the  interfacial 
zone  is  taken  into  account.  At  each  point  of  the  interface,  the  material  response  in  both  the 
normal  and  tangential  directions  is  considered  to  be  that  of  a  standard  linear  solid  (SLS). 
The  SLS  qualitatively  represents  the  behavior  of  an  idealized  cross-linked  polymer  and  can 
be  viewed  as  a  spring  in  parallel  with  a  spring  and  a  dashpot.  The  normal  and  tangential 
tractions  Tn  and  Tt  are  then  related  to  the  relative  displacement  differences  [un]/  and  [«,]/ 
across  the  interface  by 

7i  +  ^  =  (k/ +kl2)\ul]  +—(«,!  (13) 

t  r 

where  the  subscript  i  is  replaced  by  n  and  t  respectively.  As  for  the  linearly  elastic 
interface,  the  normal  and  tangential  stiffness  parameters  are  taken  to  be  equal.  Thus, 

kj  =  k/  and  k„  =  k? . 

For  convenience,  the  time  constant,  T,  in  expression  (13)  is  taken  to  be  unity  throughout. 


3.  FINITE  ELEMENT  IMPLEMENTATION 

In  order  to  solve  the  now  formulated  boundary  value  problem,  the  finite  element 
method  is  employed.  Throughout,  we  assume  small  displacements.  The  strain 
displacement  relation  then  takes  the  usual  form 

£y  =  (Uij  +  UjJ/2  (14) 

and  the  equilibrium  equation  is  written  as 

Oiy  ,j  =  0.  (15) 

For  the  fiber-reinforced  composite  with  an  interfacial  zone,  the  Principle  of  Virtual 
Work  is  written  as 

Vijfeij  dQ  +  8(p  dS  =  jf-Ti8ui  dT  (16) 
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where  Q denotes  ihe  interior  of  the  trapezoidal  region  shown  in  figure  2,  /  'is  the  external 
traction  boundary,  and  S  is  the  interfacial  traction  boundary.  The  8ui  are  the  kinematically 
admissible  displacements  (satisfying  the  periodic  displacement  boundary  conditions 
outlined  in  section  2  and  vanishing  on  the  prescribed  displacement  boundary).  The  second 
term  in  (16)  is  the  virtual  work  of  separation  of  the  matrix  and  fiber  phases,  i.e. 

5(f>  =  TnS[un\i+T;8[ut}I.  (17) 

Note  that  for  the  case  of  a  linearly  viscoelastic  interface,  the  time  dependence  of  the 
interface  enters  into  the  formulation,  and  (16)  must  be  discretized  in  both  space  and  time. 
A  detailed  formulation  of  this  discretization  is  found  in  Appendix  B. 

In  order  to  approximately  satisfy  the  impenetrability  constraint  (12)  full  Newton 
Raphson  equilibrium  iteration  is  employed  with  a  penalty-like  stress  update  scheme  in 
which  the  interface  is  taken  to  have  a  suitably  high  normal  stiffness  parameter  in 
compression. 


4.  MACROSCOPIC  BEHAVIOR 

The  numerical  results  of  the  present  analysis  provide  a  direct  method  for  calculating 
the  effective  transverse  properties  of  the  fiber-reinforced  composite.  For  the  composite 
with  a  linearly  elastic  interface,  the  numerical  results  for  the  effective  transverse  bulk  and 
shear  moduli  tire  compared  with  the  results  of  Hashin  (1989)  where  the  composite  cylinder 
assemblage  (CCA)  model  is  employed  to  obtain  the  effective  transverse  bulk  modulus,  and 
the  generalized  self  consistent  scheme  (GSCS)  model  is  used  to  determine  the  effective 
transverse  shear  modulus.  It  is  noted  that  the  CCA  model  has  been  introduced  in  Hashin 
and  Rosen  (1964)  and  generalized  to  transversely  isotropic  fibers  and  matrix  in  Hashin 
(1979).  In  the  context  of  perfect  interface  conditions  the  GSCS  model  has  been  applied  by 
Christensen  and  Lo  (1979)  to  obtain  the  transverse  shear  modulus.  For  the  linearly 
viscoelastic  interface,  the  numerical  results  are  compared  with  approximate  expressions 
derived  analytically  for  the  time  dependent  transverse  properties  of  the  composite. 

4. 1  Linearly  Elastic  Fiber/Matrix  Interface 

It  can  be  shown  analytically,  that  due  to  its  inherent  elastic  symmetry,  the  hexagonal 
array  composite  with  linearly  elastic  constituents  is  transversely  isotropic,  see  Love  ( 1927), 
Lekhnitskii  ( 1963).  'Therefore,  we  assume  that  the  effective  elastic  moduli  C* tJid  ate  related 
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to  the  Cartesian  components  of  the  macroscopic  or  average  stress  and  strain  tensors  ot] 
and  £,j  in  the  form 


Oi,  =  C 


*  — 


ijkl  £  kl. 


08) 


and  that  this  relation  takes  the  form  (7)  where  the  quantities  in  each  expression  are  replaced 
by  their  respective  averaged  or  macroscopic  counterparts  (eg.,  £ / /  is  substituted  for  £.'/ / , 
K[*  for  K j  etc.).  We  note  that  when  the  impenetrability  constraint  (12)  is  imposed,  the 
response  of  the  composite  is  observed  to  deviate  from  transverse  isotropy  when  the 
interface  stiffness  is  low,  particularly  when  it  is  compared  for  tensile  and  compressive 
loading.  For  purely  tensile  loading  conditions,  this  deviation  is  slight  however. 


Numerical  Procedure.  Under  the  condition  of  CPD  loading,  the  effective  transverse 
shear  and  bulk  moduli  of  the  composite  are  obtained  in  the  following  manner.  First,  a 
constant  displacement,  is  applied  along  the  external  boundary  EF  and  the  displacement, 
-A2,  is  applied  along  the  external  boundary  AB  (see  fig.  2).  The  solution  of  the  finite 
element  equations  yields  the  unknown,  constant  displacement  Aj  along  the  external 
boundary'  AF.  Die  average  strain  in  the  .ty  direction  is  then  given  by 
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The  average  strain  in  the  X2  direction  is  given  by 


-4  A2 
822  =  Vi  b 


(20) 


By  employing  the  relation  (5),  the  finite  element  solution  also  yields  the  remote  applied 
stress.  This  is  equal  to  the  average  stress  component  On  of  the  composite.  Since  there  is 

no  applied  loading  in  the  *y  direction,  the  average  stress  component  5y  /  -  0.  Finally,  by 
substituting  these  averaged  quantities  into  the  plane  strain  constitutive  relation  (7),  the 
effective  transverse  shear  and  bulk  moduli  are  found  to  be 
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Effective  l  ransver.se  Bulk  Modulus.  The  numerical  results  for  the  effective 
transverse  bulk  modulus,  K /* ,  are  compared  with  the  results  of  Ilashin  (19X9)  where  the 
CCA  model  is  employed.  The  CCA  expression  is  given  as 


Ej*  =  Kj, 


+  — 


Klf  -  K T 


V, 


Vr 


K 


G  / 


where  K;m  and  Gjm  are  the  transverse  bulk  and  shear  moduli  of  the  matrix,  aikH’,„  and  \  ■ 
are  the  matrix  and  fiber  volume  fractions.  The  quantity,  K;j ,  is  the  equivalent  transverse 
bulk  modulus  of  the  fiber- interlace  combination  and  can  be  expressed  as 


K  j*K_rj_kn_ 
,J  '  akn  +  2  Kjf 


(2  1) 


where  kn  is  the  normal  interlacial  stiffness  parameter,  Krt  is  the  transverse  bulk  modulus 
of  the  fiber,  and.  a,  is  the  fiber  radius.  The  numerical  and  CCA  results  for  the  effective 
transverse  bulk  modulus  are  shown  in  figure  3.  The  effective  modulus  K;  *  is  plotted 
versus  normalized  interfacial  stiffness.  The  normalization  is  chosen  such  that 


where  a  is  the  radius  of  the  fiber.  As  expected,  when  the  impenetrability  constraint  (12)  is 
not  enforced,  the  numerical  and  CCA  results  virtually  coincide  over  the  entire  spectrum  of 
interface  stiffness  parameters.  However,  when  this  constraint  is  enforced,  the  numerical 
results  deviate  from  the  CCA  results  at  relatively  low  interfacial  stiffness  parameters.  The 
magnitude  of  this  deviation  depends  on  the  loading  condition  considered  in  the  numerical 
procedure.  It  is  recalled  that  in  the  present  analysis  Kj*  is  determined  by  subjecting  the 
hexagonal  array  model  to  tensile  uniaxial  loading.  It  is  noted  that  when  the  condition  (12) 
is  not  imposed,  the  numeiicai  value  obtained  for  K p*  is  independent  of  the  loading 
condition  considered. 
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;iii  l'run>  \ .  >  vr  vV<j.*'  Modulus  As  pointed  out  hv  I  fashin  (  1 989  g  the  eftec  to  c 
transverse  shear  modulus  cannot  m  general  be  determined  by  employing  the  equivalent 
fiber -interface  concept.  ' Alike  the  CCA  expression  for  the  effective  transverse  bulk 
modulus  which  is  completely  independent  of  the  tangential  interface  stillness  parameter,  the 
effective  tiansverse  shear  modulus  is  highly  dependent  on  both  stiffness  components. 
However,  for  the  case  when  both  the  normal  and  tangential  stiffness  parameters  are  taken 
to  be  e,.u.b,  it  is  found  that  the  expression  for  the  effective  transverse  shear  module 
obtained  hv  subsuming  approximate  models  tor  the  fiber-interface  shear  modm..  ,> 
available  pence;  interface  expressions  selves  as  a  good  approximation  over  a  large  rang'-  <  •: 
mte: laciiti  'tit!  less  values. 

m  the  present  analysis,  we  employ  an  approximate  expression  tor  the  shear  momma- 
of  the  fi'v  m-erface  combination  which  is  written  as 

uC’/gf, 

v.  /if  -  ~  r  "  %  '  v  •  .  1 

(l  fci  4-  J  (  ’  i 


where  k,  l*  the  tangential  inter! acini  stiffness  parameter,  (•:■  is  the  trait*  verse  shear 
module-,  of  the  fiber  phase,  and  «,  is  die  fiber  radius.  It  is  noted  that  when  the  normal  and 
tangential  stiffness  parameters  are  taken  to  be  equal,  the  result  obtained  by  substituting  ;2n, 
Vi  the  transverse  shear  modulus  of  the  fiber  in  the  expression  ot  Christensen  and  l.o 
i  !97‘>;  is  indistinguishable  from  the  CISC'S  result  of  Hashin  ( 1989).  In  figure  -C  the  ( iSCS 
results  obtained  by  Hashin  (19.SU)  are  compared  with  the  numerical  results.  Results 
o.-nau'ed  hv  substituting  Cm  for  (he  shear  modulus  of  the  fiber  phase  in  the  lower  bound 
expression  fo:  arbitral y  phase  geometry  given  by  Hashin  (1979)  are  also  shown  in  the 
figure  The  effective  transverse  shear  modulus  normalized  with  respect  to  the  shear 
modulus  of  the  matrix  is  plotted  versus  normalized  mterfaeial  stiffness.  Again,  the 
normalization  C  chosen  as  m  (25).  When  the  impenetrability  constraint  (12)  is  enforced, 
the  numeric  ai  results  deviate  significantly  from  the  GSCS  results  at  relatively  low  intcrfacia! 
stiffness  parameters.  When  the  constraint  ( 12)  is  not  enforced,  the  deviation  at  low 
stif tncs  >cs  is  less  significant.  Again,  when  ( i2)  is  enforced,  tfic  magnitude  oi  the  deviation 
depends  on  the  loading  condition  imposed  on  the  hexagonal  array  model  in  the  numeric.il 
proc-v  It  o  noticed  that  under  tliese  circumstances  the  modified  lower  bound 
expression  serves  as  a  good  model  for  the  effective  transverse  shear  modulus  of  the 
composite.  The  corresponding  results  for  the  effective  transvei.se  bulk  and  she..*  modulus 
versus  f  hep  volume  fraction  for  a  fixed  normalized  interfacial  stiffness  parameter  o n  j ) 
am  -mown  in  figure  V 
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4.2  l., nearly  \  Yu .  t  laste  ribcr  Matrix  Interface 

In  principle.  n  expressions  for  the  effective  moduli  of  the  composite  with  a  lincasl> 
elastic  interlace  ere  known,  the  correspondence  principle  can  be  employed  in  order  to 
calculate  the  effective  time  dependent  expressions  for  the  composite  with  a  linearis 
viscoelastic  interface.  1  nfortunately,  the  exact  elastic  expression  for  the  composite  under 
consideration  remains  to  be  determined.  However,  it  is  noticed  that  the  CCA  results  for  the 
effective  transverse  bulk  modulus  given  in  the  previous  section  compare  closely  with  the 
numerical  results  for  a  wide  range  of  stiffness  parameters  over  the  entire  range  of  volume 
fraction.*  considered  We  thus  choose  the  linearly  elastic  CCA  expression  and  employ  the 
correspondence  principle  to  obtain  a  model  for  the  time  dependent  effective  modulus. 
A' /*(?'.  of  tKc  fiber  reinforced  composite  with  a  linearly  viscoelastic  interl'acia!  /one.  In 
order  to  obtain  an  analytic  expression  for  the  time  dependent  effective  tiai'svcoe  he.!: 
modulus.  G  ■  ’it),  the  mod  tried  lower  bound  expression  obtained  in  the  previous  -cctio.i 
chosen  since  it  is  algebraically  simple  and  lends  itself  well  to  analytic  Laplace  ‘larsfo:;:! 
inversion.  However,  it  is  noted  that  the  usefulness  of  the  resulting  time  dependent  m/x.Y: 
for  the  effective  transverse  shear  modulus  is  limited  since  it  is  only  good  for  the  special 
case  when  the  normal  and  tangential  interfacial  stiffness  parameters  are  assumed  to  he 
synchronous  during  interfacial  stress  relaxation. 


Analytical  Model.  By  substituting  the  transformed  relaxation  function  of  the  SI  S. 
.h7(,y),  into  (24)  and  (26),  the  transformed  expressions  sKf/U)  and  .vG/yriv)  are  obtained. 
The  substitution  of  these  expressions  (orKjj  andG/y  in  the  CCA  result  (23)  and  the  lower 
bound  result  of  Kashin  (1979)  yields  sK;*(s)  and  sGj*(s).  Generally  a  problem  arises 
here  in  inverting  back  into  the  time  domain.  However,  for  this  relatively  simple 
viscoelastic  interface  model,  the  relations  for  Kf*(s)  and  Gj*(s)  are  easily  inverted  yielding 
a  model  for  the  time  dependent  effective  moduli  Kj*(t)  and  Gj*{t).  The  effective  relaxation 
function  in  bulk  can  be  expressed  as 

Krit)  =  K rJ  +  (K,/  -  KtJ  )e~tllK  > 


while  the  effective  relaxation  function  in  shear  can  be  written  as 


Grit)  =  (- 


(Grg*  -  Gj„*  )e 


l!Tr. 


\  f  X  i 


where  K  ru*  and  G'/C  are  the  glassy  or  instantaneous  values  for  the  effective  n 
shear  and  bulk  moduli  of  the  composite,  K  c  *  and  Gj.-*  are  the  corresponding  ! 


values  obtained  as  time  approaches  infinity,  and  TK  and  Ta  are  the  respective  time 
constants  in  Lansverse  bulk  and  shear.  These  quantities  are  defined  in  Appendix  A. 


Comparison  of  Analytical  Model  With  Numerical  Results.  The  numerical  procedure 
employed  to  determine  the  effective  relaxation  moduli  Kj*(t)  and  Gf*(t)  is  almost  identical 
to  that  used  in  the  numerical  determination  of  the  elastic  transverse  moduli.  The  step 
displacement  A2H  (t)  is  applied  along  the  external  boundary  EF,  and  -zD//(f)  is  applied 
along  tire  external  boundary  AB  (see  fig. 2).  Using  the  relations  (21)  and  (22),  the  effective 
moduli  are  computed  and  plotted  as  a  function  of  time.  The  effective  relaxation  moduli 
Kj*(t)  and  G  f  it)  predicted  by  the  analytical  model  and  those  obtained  numerically  are 
plotted  for  a  fixed  fiber  volume  fraction  ( Vj  -0.5 )  as  shown  in  figure  6.  The  normalized 

l  ■> 

stiffness  parameters  k  and  k~  are  varied  such  that  the  glassy  stiffness  of  the  interface 
ranges  from  a  normalized  value  of  ( kg=10 )  to  a  normalized  value  of  (kg-J).  Throughout, 
the  ratio  of  glassy  to  long-term  interfacial  stiffness  is  taken  to  be  (kg/k^=  10).  The 
normaliz.ation  is  chosen  such  that 


(29) 


where  a  is  the  fiber  radius  and  the  superscript  i  takes  on  the  values  l  and  2  referring  to  the 
respective  chains  of  the  SLS.  As  shown  in  the  figure,  the  approximate  analytical  model 
compares  extremely  well  with  the  numerical  results  over  the  range  where  the  interfacial 
stiffness  is  relatively  high.  As  expected,  when  the  constraint  (12)  is  enforced,  for  the 
relatively  low  glassy  stiffness  value  ( kg=l ),  the  numerical  results  lie  below  the  derived 
results  for  the  transverse  bulk  relaxation  modulus,  and  above  the  derived  results  for  the 
transverse  shear  relaxation  modulus.  It  is  interesting  to  note  that  when  the  impenetrability 
constraint  (12)  is  relaxed,  the  numerical  results  virtually  coincide  with  the  model  over  the 
entire  range.  These  results  for  the  Mid-CPD  loading  direction  do  not  differ  significantly. 


5.  MICROSCOPIC  BEHAVIOR 

In  this  section,  we  turn  from  a  macroscopic  perspective  and  focus  on  the  local  or 
microscopic  behavior  of  the  transversely  loaded  fiber-reinforced  composite.  In  particular, 
we  examine  the  influence  of  interfacial  stress  relaxation  on  the  stress  fields  in  the  matrix 
material  contiguous  to  the  interface.  Throughout,  the  composite  is  subjected  to  a  step 
displacement  in  the  Mid-CPD  loading  direction,  the  more  severe  of  the  two  loading  cases 
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considered,  and  normalized  circumferential  stress  distributions  in  the  matrix  material  just 
outside  of  the  interface  are  obtained  during  the  numerical  simulation  of  stress  relaxation. 

During  the  interfacial  stress  relaxation,  the  circumferential  stress,  oefh  normalized  with 
respect  to  the  remote  applied  stress,  a »  is  plotted  versus  angle,  6,  as  shown  m  figure  7. 
In  the  upper  diagram,  the  normalized  glassy  stiffness  of  the  composite  is  taken  to  be 
( kg=I0 )  and  the  ration  of  glassy  to  long  -term  stiffness  is  taken  to  be  (kg/k^=10).  Again, 
the  time  constant,  t,  is  taken  to  be  unity  for  convenience.  As  shown  in  the  figure,  the 
normalized  circumferential  stress  increases  substantially  as  time  progresses  and  as  the 
effective  stiffness  of  the  interface  decreases.  Note  that  as  the  interface  relaxes,  a  relatively 
abrupt  change  in  the  distributions  occurs  at  approximately  75  degrees.  This  is  due  to  a 
local  stress  concentration  at  the  transitional  point  where  the  compressive  region  of  the 
interface  begins.  The  resulting  distributions,  as  interfacial  stress  relaxation  proceeds,  are 
shown  in  the  lower  diagram  where  the  normalized  glassy  stiffness  of  the  interface  is  taken 
to  be  (kg=l).  As  time  progresses,  a  noticeable  shift  occurs  in  the  distributions.  The 
maximum  normalized  circumferential  stress  gradually  relocates  from  approximately  75 
degrees  to  approximately  the  45  degree  mark.  In  a  recent  experimental  study  by  Daniel  et 
al.y  (1989),  when  a  silicon  carbide/glass-ceramic  composite  having  closely  spaced,  near 
hexagonally  packed  fibers  was  transversely  loaded,  fracture  was  observed  to  initiate  in  the 
form  of  radial  cracxs  in  the  matrix  at  this  45  degree  location. 

6.  CONCLUDING  REMARKS 

In  the  present  work,  the  transverse  loading  of  the  hexagonal  array  composite  is 
examined  firm  both  a  macroscopic  and  a  microscopic  point  of  view.  In  the  macroscopic 
analysis,  the  composite  is  subjected  to  tensile  loading  and  the  effective  transverse  properties 
are  obtained.  For  the  case  of  a  linearly  elastic  interface,  the  numerical  results  are  compared 
with  expressions  obtained  from  the  composite  cylinder  assemblage  and  generalized  self 
consistent  scheme  models.  It  is  found  that  when  the  impenetrability  constraint  is  imposed, 
the  mechanical  behavior  of  the  hexagonal  array  composite  deviates  slightly  from  transverse 
isotropy  and  the  numerical  results  deviate  significantly  from  the  analytic  expressions  at  low 
interfacial  stiffness  values.  The  magnitude  of  this  deviation  depends  on  the  loading 
configuration  considered  in  the  numerical  procedure.  Note  that  we  impose  this  constraint 
to  avoid  an  unrealistic  interpenetration  of  the  matrix  phase  into  the  fiber  phase.  For  the  case 
of  a  linearly  viscoelastic  interface  the  numerical  results  for  the  time  dependent  transverse 
moduli  are  compared  with  approximate  expressions  derived  analytically.  Good  agreement 
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is  found  among  the  results  over  a  large  range  of  interfacial  stiffnesses.  In  the 
micTomechanical  analysis,  normalized  circumferential  stress  distributions  are  obtained  in 
the  matrix  material  contiguous  to  the  interface  during  interfacial  stress  relaxation.  It  is 
found  that  these  distributions  change  substantially  as  the  effective  stiffness  of  the  interface 
decreases.  The  maximum  circumferential  stress  concentration  is  found  to  occur  at 
approximately  45  degrees  at  relatively  low  interfacial  stiffness  values. 
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APPENDIX  A:  DEFINITIONS 


In  this  appendix,  we  define  the  quantities  contained  within  the  derived  expressions 
(27  and  28)  for  the  time  dependent  relaxation  moduli.  The  expressions  for  Kjg*  and  Gjg* 
are  obtained  by  making  the  following  substitutions: 

I  2 

1.  Substitute  for  kn  in  (24)  the  sum  ( kn  +  kn  )  of  the  stiffness  parameters  of  the  SLS 

obtaining  the  intermediate  expression  Kjg  for  the  glassy  effective  bulk  modulus  of  the 

1  1 

fiber-interface  combination.  Similarly,  substitute  for  kt  in  (26)  the  sum  (k,  +  kt~ )  for 
the  glassy  effective  shear  modulus  Gfg- 

2.  Substitute  the  expression  Kjg  into  the  CCA  result  (23)  and  Grg  into  the  lower  bound 
expression  for  arbitrary  phase  geometry  of  Hashin  (1979). 


Thus,  the  effective  glassy  moduli  are  written  as 


Kjg*  =  Kfm  + 


GtS*=  Gfm  + 


Vf 


Kjg  Kt  m 


Krm  +  Gjr, 


il 


Gto  -  Gt, 


{Kjm  +  2Grm)V m 
2GTm  ( KJm  +  GIm). 


(Al) 

(A 2) 


where 

(A3) 

(A4) 

and  the  long-term  effective  moduli  are  expressed  as 

Kr~*  =  KTm  +  - 7 - ^ - v -  (A5) 

_ i _  +  _ v  m _ 

Kt co  -  KTm  Klm  +  Gjm 


l  2 

p  -  aKTf(kn  +  kn  ) 
^Tq  —  12 

a(kn+kn)  +  2KTf 


r  _  aGTf(kt V  kt2) 

\JTg  —  12 

a  (k,  +  k,  )+ 2GTf 
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GjJ=  Gj„ 


1 


Gt°°  ~  Gjm 


+ 


L _ 

(Kjm  +  2GTm)Vm 

2  Gjm  C^Tm  +  Gfm) 


(A6) 


where 

kt~  =  -aKrthL  (A7) 

akn‘  +  2KTf 

Gr«=  -afTfk‘ — .  (A8) 

akt  +  2Gjf 


Finally,  the  time  constants  TK  and  TG  are  written  in  terms  of  the  time  constant,  r  ,  of  the 
SLS  as 


J  (kr!  +  kn)  1+  m  -  Kim) 

kn  KTg  1+  m  (KTg-  KTm) 

j  (&/  +  r^°°  1+  n  (Gt~  ~  GTm) 

kt  Grg  1+  n  ( Gjg  -  Gjm) 

where 

*Tm  +  Gfm 


(A9) 


(A  10) 


(All) 


(KTm+  2GTm)Vm 
2Gjm  (Kt m+  Grm )  ' 


(A  12) 


APPENDIX  B:  FINITE  ELEMENT  FORMULATION  FOR  LINEARLY 
VISCOELASTIC  INTERFACE 

For  the  case  of  a  linearly  viscoelastic  interface,  the  time  dependence  of  the  interface 
enters  into  the  formulation,  and  the  Principle  of  Virtual  Work 

Oijdeij  dQ  +  J5  8<p  dS  =  jrTi8ui  dr  (Bl) 
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must  be  discretized  both  in  space  and  time.  The  following  implementational  details  are 
based  on  a  solution  procedure  that  has  been  discussed  by  Taylor,  Pister  and  Goudreau 
(1970).  Writing  (Bl)  in  incremental  form  at  time  tn+i  yields 

j"^  AOijSe^  d.C2  +  ATi8[ui\i  dS  = 

jpATM  dr-  [  (TijSEij  dQ  +  J5  TiS[Ui],  dS  -  jpTrfut  dr  ]  „  (B2) 

where  (<7y)n+;  =  (<Ty)„  +  (4cry)„  and  (T,)„+;  =  (7/)B  +  (ATt)n.  Using  (B2)  as  our  starting 
point,  we  now  express  the  components  7,  of  the  interfacial  traction  vector  in  terms  of  the 
[«,]/  as  follows.  The  standard  linear  solid  can  be  viewed  as  being  comprised  of  two 
Maxwell  elements  in  parallel,  one  of  which  has  an  infinite  time  constant  r.  The  governing 
differential  equation  for  one  Maxwell  element  is  simply  expressed  as 

T  i 

T/+  =  (£/)[«;]/  (B3) 

xJ 

where  the  superscript  j  refers  to  the y,/,  Maxwell  element  of  the  SLS.  For  extension  at  a 
constant  rate,  [tij/,  the  exact  solution  of  (B3)  is 

T/(t )  =  (k/mhr^l-exp  (-t  lx)}.  (B4) 

Letting  ( A[ui]i)n  =  ([«;]/)„+/  -  ([«,]/)„  and  At  =  tn+j  -  tn,  the  interfacial  traction  components 
7,  at  time  tn+j  are  obtained  by  expressing  (B4)  in  incremental  form  and  summing  the  partial 
stresses.  Thus 

2  r  -i 

(7«  )n+i  =  I  L  T-exp  (- At  /  x)  +  k/ A[ui\i  h(At)  J  „  (B5) 

j=l 


where 

h{At)  =  x \l-  exp  {-At  /  t;)  )  /  At .  (B6) 

Note  that  as  the  time  constant,  xJ,  approaches  infinity,  h(At)  in  (B6)  approaches  unity. 
Referring  back  to  equation  (B2),  we  introduce  standard  finite  element  shape  functions,  N, 
and  use  the  arbitrariness  of  the  variations  in  the  nodal  displacements  to  obtain  the  discrete 
form 
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(B7) 


S  [  La/dA  ^  +  Js  H*  dS  ]  „+;  = 

e=/  J 

I  [  \rNeTTn+ldr-  \SBrTn+1dS-  [QBTGndQ  ] 
e-1  J  *  j 

where,  the  vector  Ade  contains  the  discrete  nodal  displacement  increments  and  ne  refers  to 
the  number  of  elements  in  the  model.  The  first  term  on  the  left  side  of  (B7)  contains  the 
symmetric  gradient  operator  matrix  Be  and  the  plane  strain  constitutive  matrix  De.  The 
second  term  on  the  left  side  of  (B7)  constitutes  the  interface  stiffness  contribution  of  the 
interface  layer,  see  Gosz  (1989)  for  interface  element  formulation.  Here,  the  interface 
constitutive  matrix  De  is  a  2  x  2  diagonal  matrix  whose  components  are  given  by 

2 

Du  =  I  h(At)knJ  (B8) 

j=l 

2 

D22  =  I  h(At)k/ 

hi 

where  kj  and  k/  are  the  normal  and  tangential  stiffness  parameters  associated  with  th cjth 
element  of  the  SLS  and  h(At)  is  defined  by  (B 6).  Note  that  this  entire  procedure  can  be 
easily  extended  for  more  general  constitutive  relations,  for  example,  by  modeling  the 
normal  and  tangential  material  response  of  the  interface  with  a  greater  number  of  Maxwell 
elements. 
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Table  1.  Flastic  constants  of  the  graphite/epoxy  material  system  considered. 


EA(Gpa  ) 

Ej  ( Gpa ) 

Vr 

GA(Gpa ) 

Gj  ( Gpa ) 

Kj  ( Gpa  ) 

Graphite  fiber 

232 

0.279 

15.0 

0.490 

24.0 

5.03 

15.0 

Epoxy  matrix 

5.35 

0.354 

5.35 

0.354 

1.976 

1.976 

6.76 

Fig.  1.  Mid-CPD  loading  configuration  for  the  hexagonal  array  composite. 


Fig.  2.  Schematic  of  the  representative  cell  of  the  hexagonal  array  composite. 
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Fig.  4.  Normalized  effective  transverse  shear  modulus  versus  interface  stiffness. 
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Fig.  6.  Comparison  of  relaxation  model  in  transverse  bulk  and  shear  with  the  numerical  results 

(linearly  viscoelastic  interfacial  zone). 


